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ABSTRACT 

I describe how the states of a discrete automata with p sites, each of which may 
be off or on, can be represented as Majorana spinors associated to a spacetime with 
signature (p,p). Some ideas about the quantization of such systems are discussed 
and the relationship to some unconventional formulations and generalizations of 
quantum mechanics, particularly Jordan's spinorial quantum mechanics are pointed 
out. A connection is made to the problem of time and the complex numbers in 
quantum gravity. 
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Introduction 

Andre Trautman has had, throughout his career, the happy knack of realizing long 
ahead of others what will ultimately turn out to be just the right language and 
concepts needed to most elegantly and economically encapsulate our emerging ideas 
about the physical world. Such languages must be both precise and universal, easily 
understood and yet sufficiently subtle as to suggest new avenues for exploration 
and exploitation. Often called geometrical, they may in fact incorporate profound 
algebraic insights as well. 

I first became aware of Andre's talent for clear structured exposition of a physi- 
cal theory on reading as a student his 1964 Brandeis Lectures, appropriately entitled 
" Foundations and Current Problems of General Relativity". A later notable ex- 
ample of his prescience was his insight that fibre-bundle theory would eventually 
prove to be an everyday tool for physicists. More recently he has been exploring 
the world of spinors and geometrical algebra. In this article I wish, among other 
things, to indicate how this might just turn out to be the appropriate language to 
describe such disparate subjects as the spread of forest fires, neural nets and the 
problem of memory or the development of avalanches, just as it has already shown 
its worth in statistical mechanics. 

My own interest was first aroused by some ideas of 't Hooft's on cellular au- 
tomata and quantum mechanics [1,2]. During a stay in Utrecht I tried to understand 
what was the connection bwetween "spin models" and spinors which was so effec- 
tively exploited by by Onsager 50 years ago in his solution of the Ising model. Since 
then I have been struck by how often the same ideas seem to crop up in many other 
areas of physics of which I have only a superficial knowledge. Strange analogues of 
quantum mechanics, quantum field theory and the Dirac equation seem to arise all 
over the place. Clearly there is something universal going on here. What follows 
is my attempt to understand what that might be. I do not claim that I am alone 
in realizing this connection (Dubois- Violette, for example, has described similar 
things [3]) but my emphasis on the essential reality (in the mathematical sense) of 
the construction and my attempts (in the final section) to relate the discussion to 
some ideas in quantum gravity is perhaps new. In any event I dedicate it to Andre 
hoping that he will find my exposition of interest and that it may perhaps provoke 
him to look further if he hasn't already done so. 
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1. Master Equations and Partition Functions 

The common theme of many of the applications I alluded to is a Markov process 
in which the time development of a "probablity vector" p = (pi,p2, ■ ■ ■ ,p n ) evolves 
according to a " Master Equation" of the form 

= {~ w ijPi + w HPi) ( L1 ) 
= J2 w ijPv (1-2) 

3 

where the "transition rates" iOy are positive and constrained by the requirement 
that 

so that for example 

Wll = -W12 - W13 • • • - Win- (1-4) 

Geometrically speaking (1.3) tells us that the probability vector lies on the (n — 1)- 
simplex E n _i in the positive orthant WJ of M n * satisfying 

i.p = l (1.5) 

where i = (1,1, .. .1) is the un-normalized completely ignorant probability vector 
obtained by taking all a priori probablities to be equal. The set of probability vec- 
tors, also called states, is convex. The extreme points of the state space £ n _i, may 
be are called "pure states", and correspond to the vertices of the simplex, that is 
to the unit basis vectors e; in IR n and the remaining, non-extreme states are called 
" mixed" . Note that the language of pure and mixed states is as applicable to the 
classical systems being considered here as it is to quantum mechanical systems. In- 
deed according to one view point what makes the present set up classical is precisely 
the fact that the space of mixed states M. is a simplex [4,5,6] and the space of pure 
states V are its vertices. 

* By M n , I mean the standard n-dimensional vector space equipped with its 
positive definite inner product. I shall not, therefore, distinguish between vectors 
and covectors. All indices will therefore be lowered. 
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The solution of (1.1) may be written as as 

Pi(t) = exp(twij)pj(0) (1.6) 

= $yPi(0), (1.7) 

where $ = exptw is a stochastic matrix, i.e a square matrix all of whose rows sum 
to unity. 

The stationary, i.e time-independent, states, or equilibria have probabilities 
p equ ■ In thermodynamic applications each pure state e$ is assigned an energy Ei 
one has at temperature T = 

P ^ = eM-m/m (i.8) 

where: 

Z(/3)=^exp(-^). (1.9) 

i 

is the partition function. Given the energies Ei the challenge is to calculate the 
partition function Z{(3). This was achieved by Onsager [7] using Pauli matrices 
including the pure iamaginary matrix oi for the two-dimnsional lattice Ising model. 
I shall argue shortly that, despite appearances and as befits the purely classical 
nature of the problem, these spinorial techniques actually involve the real, i.e. Ma- 
jorana spinors of a Kleinian spacetime with as many time coordinates as it has space 
coordinates. 

More generally given an observable O which takes values O n , its expectation 
value (O) at time t is given by 

(o) = X>Oi. (i-io) 

i 

One may think of the classical (commuting) observables as matrices acting on 
R n which happen to be diagonal in the basis provided by the pure states. One may 
then re-write (1.10) as 

(0) = i.O P . (1.11) 
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2 Quantizing Discrete Automata 

I have already commented upon the superfical resemblance of this classical 
formalism to standard quantum mechanics. If we really were dealing with quantum 
mechanics, then the corresponding configuration space Q would be just the the 
discrete set of n points Q n = {e^} and the Hilbert space H would be L 2 (Q n ,C) 
which may think of as the standard n-dimensional Hermitean vector space C n = 
R n <g> C with orthonomal basis given by {ej}. 

Thus even though the (imear) equation (1.2) is like the (imear) Schrodinger 
equation we face the difficulty that 

(i) The probablity vector p has components which are real and positive. 

(ii) The probablity vector p is normalized using the L 1 norm (1.4) rather than the 
usual L 2 norm, 

One might consider, as does Wheeler in considering thesis topics [9] , passing to real 
quantum mechanics [8] by taking a square root, i.e. introducing the vector 

V> = ((Pl)*,(P2)V..(Pn)*)- (2- 1 ) 

The normalization condition (1.4) now tells us that ip lies on the unit n-sphere 
S n ~ 1 C M n . However now the equation of motion (1.2) becomes non-linear. More- 
over there is some ambiguity in taking the square roots in (2.1). If we take them all 
to be positive we map the simplex £ n _i onto the positive orthant of the (n — 1)- 
sphere or projectively speaking, the interior of an (n — l)-simplex sitting inside real 
projective space RP n_1 . In real quantum mechanics one usually thinks of the pure 
states as all real vectors ip except that one identifies those which differ by a non-zero 
real multiple. Thus the the pure states roam all over IRP n_1 . To restrict them to 
a part of IRP n_1 would in effect put a restriction on the Superposition Principle : 
only postive real combinations would be allowed. 

The idea pursued by 't Hooft is different. He considers 

(i) The evolution to be reversible and time to be discrete 

(ii) Each site to be occupied with certainty. 

Thus we have a classical reversible discrete automaton whose time evolution is given 
by iterating a permutaion matrix U : Q n — > Q n which permutes the basis elements 
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{ej of W 1 . After N steps we obtain U N . Clearly U N is a special case of an 
orthogonal matrix 

U N cSO(n;R) (2.2) 

or indeed a unitary matrix 

U N cU(n;C). (2.3) 
't Hooft's idea is to find a Hermitean matrix H such that 

U N = exp(-iNH). (2.4) 

The Hermitean matrix H is to be thought of as the Hamiltonian of the quantum 
mechanical system associated to the original classical automaton. 

Quite aside from the ambiguity in taking the logarithm in order to obtain H, 
there are two rather puzzling things about this. 

(i) Usually we only consider classical first order equations governing the motion in 
classical state space or phase space V which is even dimensional. The motion 
in configuration space Q is usual governed by second order equations. Thus we 
might have tried to identify V with the discrete set of states {e^}. 

(ii) Usually the quantum mechanical Hilbert space Ti is not given by L 2 (V), rather 
one seeks a prescription (e.g a "polarization") for cutting L 2 (V) down to a 
Hilbert space of functions depending upon half as many variables. The standard 
continuous case is of course when the phase space is the tangent bundle of the 
configuration space V = T*(Q). 

This suggests to me that one should restrict attention to the case when the 
number of discrete states n is even. One could then seek to embed the time evolution 
operator in C), 

71 

U N CU(-;C). (2.5) 

In other words one should try to endow the original real vector space R n with a 
complex structure J thus allowing us to regard it as C% . Acting on the states one 
would have: 

e l -> J(eO (2.6) 

such that J 2 = — 1 and moreover one would require that this action commutes with 
the time evolution: 

UJ-JU = 0. (2.7) 
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The complex structure J would then play the analogous role to that of a " polariza- 
tion" in the geometrical quantiztion of a symplectic manifold. 

An interesting and different example where one passes from a classical to a 
quantum desription of a discrete system arises in the the theory of quantum com- 
putation [11]. Classically a single bit of information is carried by two-state system. 
Quantum mechanically a qbit is carried by a two-state quantum sytem. However 
in that case it is physically more realistic to think of the classical two state system 
arising from the quantum-mechanical two-state system in the limit that all quantum 
coherence is lost rather than thinking of the quantum system as arising from the 
"quantization of the classical system " . Thus classically one has a 2 x 2 diagonal 
density matrix whose entries are and +p| say. Thus the remarks of the last 
paragraph do not apply. 

3 Site Models 

So far I have said nothing about what the probabilities pi are supposed to be for. 
Moreover the systems we have considered have no " spatial" structure and so notions 
like locality are not defined. In many applications one considers a certain number, 
p of sites which may be occupied or not. Typical examples are "classical spins" 
which may be "up" or "down" in the Ising model or trees in a forest which may 
be either "green" or "burnt out". The total number of pure states n is therfore 2 P . 
This is exactly the dimension of the Grassmann algebra A*(M P ) over MP and the 
connection arises as follows. We have 

A*(RP) = ^©AP(RP). (3.1) 

g=0 

where A P (M P ) are the (/-forms on M p . The space of (/-forms is spanned by the simple 
(/-forms. These correspond to the (/-planes through the origin of M p which contain 
q distinct basis vectors. There are q \(pL q )\ of these and they correspond to the 

ways of choosing q occupied sites. Thus for exmple A°(IR P ) = R correponds to no 
sites being occupied and A P (R P ) = M correponds to all sites being occupied. Hodge 
duality acting on A*(M P ) reverses the state of occupation of the sites. 

Thus if we think of a single site ex and p; is the probability of it being empty 
and p-f is the probability of it being full we have the probablity vector 

P = P4+PT e i- ( 3 - 2 ) 
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If two sites are involved we have 

P = Pu +Pn e i +PU e 2 +PtT e i Ae 2- (3-3) 

As an example one could take Empedocles's theory of the four elements, fire, earth 
air and water based on the two pairs of contrary attributes hot /cold and dry /wet. 
Ramon Lull's Ars Magna [10] is essentially the generalization to more than two 
sites. 

These last two examples look, on the face of it, neither quantum mechanical 
nor spinorial but let's define real creation and anhillation operators a a and a\ where 
the index a taken from the beginning of the latin alphabet run from 1 to p by their 
action on a general probability vector p e A* (W) by 

4p = e a A p (3.4) 

and 

a a P = i ea P, (3.5) 

where A is the exterior product and i denotes the interior product. It follows that 

a\a b + a b a\ = 5 ab (3.6) 



44 + a\a\ = (3.7) 

and 

a a a b + a b a a = 0. (3.8) 

In fact it is straightforward to check that, with respect to the euclidean inner product 

induced on A*(IR P ) form the euclidean inner product on R p , that a a and a£ are 
adjoints or in the natural basis transposes: 

a\ = (a a y. (3.9) 

If one changes basis by introducing 

7 ±a = a a ± aj (3.10) 

one has 

1±a1±b + l±bl±a = ±2<5 a6 (3.11) 
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and 



7 +a 7_ 6 + 7-67+a = 



(3.12) 




1} generate the 



where IR(2 P ) is the algebra of real 2 P x S p matrices. In other words we may think of 
the probability vectors p as being Majorana spinors for SO(p,p; M), hence my title. 
The space of mixed states Ai is thus the set of Majorana spinors whose components 
are non-negative and which sum to unity. The set of pure states V are those spinors 
with one non-vanishing component equal to unity. Unfortunately they do not in 
general correpsond to what are called " pure spinors" . 

Note that the Clifford group itself is the general linear group GL(2 P ;'R). The 
annihillation and creation operators {a a , a\} themselves are associated with lightlike 
directions. The volume element rj G A P (M) is given by 



Thus r\ = 1 for states with an even number of occupied sites and r\ = — 1 for states 
with an odd number of occupied sites. It follows that the direct sum decomposition: 



corresponds to the direct sum decomposition of the space of Majorana spinors into 
Weyl and anti-Weyl spinors. Weyl spinors have an even number of occupied sites 
and anti-Weyl spinors an odd number of occupied sites. 



a=p 



»7=IJ(l-2aIa a ). 



(3.13) 



a=l 




(3.15) 
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Spinors, the Light Cone and Jordan Quantum Mechanics 

In view of the connection emphasised by Penrose between 2-component Weyl spinors 
of SO(3 7 1), the lightcone of four-dimensional Minkowski spacetime and the quan- 
tum mechanics of spin It seems appropriate, while considering quantum mechanics 
and spinors, to recall a generalization of the usual quantum mechanical formal- 
ism due to Jordan which dispenses with wave functions and takes density matrices 
as being fundamental. As remarked by Townsend, [12] this is especially intrigu- 
ing because of Hawking's well known suggestion that quantum coherence may be 
lost in quantum gravity due to black hole effects. Jordan requires of his "observ- 
ables" p they satisfy the axioms of a real commutative but non-associative algebra 
A = {p, •} now called a Jordan Algebra. The irreducible finite-dimensional algebras 
were classified by Jordan, von Neuman and Wigner [13]. There are four series and 
one exceptional case: 

(i) J(M n ) 

(ii) i£ 

(iii) H% 

(iv) H* 

(v) tf 3 ° 

Of these H^, consist ofnxn real symmetric, Hermitean and quater- 

nionic Hermitean matrices respectively , corresponding to quantum mechanics over 
the fields M, C, EI for which the product • is given by 

Pi •Pi = ^(^1^2 + P2Pi)- (4-1) 

The exceptional case H® is related to 3 x 3 the octonionic matrices. The (n + 1)- 
dimensional sequence J(M n ) is based on Clifford multiplication associated to M. n . If 
{1, 7^} is an orthormal basis for the Clifford algebra Cliff (M n ) a general mixed state 
may be written as 

and the Jordan product satisfies 

li*lj = 5ij- (4.3) 
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Note that the non-associative algebra J(IR n ) is not a subalgebra of associative al- 
gebra Cliff (R n ) because they use a different multiplication law. However given a 
matrix representation of the Clifford algebra one obtains a matrix represention of 
the Jordan algebra using (4.1). 

Note that as algebra's Jordan algebras are necessarily real and indeed one may 
regard this as a part of the motivation: to deal only with real observables quantities. 
However it may happen of course that the spinorial Jordan algebras J(IR n ) based 
on the real Cliford algebra Cliff (IR) n ) are isomorphic to algebras of complex valued 
matrices. Whether or not this is true depends upon the dimension n. This happens 
for the case n = 3, which may be said in some way to account for the sucess of 
complex methods in four-dimensional general relativity. 

Let us return to the interpretation of spinorial Jordan algebras. If n; = — 

(x fe x fe )2 

then provided one sets xq = 1 every non-trivial mixed state p may be expressed 
uniquely in terms of two pure states which are associated to two primitive idempo- 
tents or projection operators E± = | (l ± n^i) 

P = pE+ + (l-p)E-. (4.4) 

The idempotents E ± satisfy 

E\ = E±. (4.5) 

and 

E±»E^ = 0. (4.6) 

Since \x\ = \p — ||, and one wishes to interpret p as a probability to be in the 
pure state associated to E +: observables are confined to lie within or on the surface 
of a ball of unit diameter in M n . The mixed states M lie in the interior and the pure 
states V, which correspond to idempotents of the algebra, lie on the unit (n — 1) 
sphere. 

XiXi = 1. (4.7) 

Geometrically speaking, a general mixed state lies on a unique diameter of the ball. 
The ends of which are the pure states. The probabilities to be in the two pure states 
are given by the distances to the endpoints of the diameter. 

Clearly we may think of the set of pure states of the Jordan algebra J(IR n ) as 
the light cone of (n + 1) -dimensional Minkowski spacetime R"' 1 . There are four 
special cases: 
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(i) J(R) = Cliff (R) = E = R © R, 

(ii) J(M 2 ) = tf 2 K , 

(iii) J(M 3 ) = tf 2 c , 

(iv) J(M 5 ) = if 2 c . 

Case (i) is that of a single a classical spin. The symbol E denotes the double 
or hypercomplex numbers. As an algebra it is reducible. The space of pure states 
is S° = Z 2 . 

Case (iii) coincides with the smallest possible non-trivial standard quantum 
mechanical system having just two states with Hilbert space C 2 . The space of pure 
states is S 2 = CP 1 . 

Case (ii) is the real quantum mechanical version of this with Hilbert space R 2 . 
The space of pure states is S 1 = MP 1 . 

Case (iv) corresponds to quarternionic quantum mechanics with Hilbert space 
H 2 . The space of pure states is S 4 = HP 1 . 

The other cases differ radically from standard quantum, or non-standard me- 
chanics. In general there is no simple correspondence beween light rays or pure 
states and spinors and the spinors are not in general Majorana. However, as we 
shall see shortly, one may also discuss master equations in the context of spinorial 
Jordan mechanics. 

It is convenient to consider only observables which have no component along 
the direction of the unit matrix. Then, up to scale, they may be identified with 
points on the sphere with unit diameter. 

If we have a matrix representation we may express the expectation value of an 
observable a in a state p as 

2 

(a) = — Trpa = 2^^, (4-7.) 

where v is the dimension of the representation of the Clifford algebra, ii.e. 2^1. We 
may write the expectation value as 

(a)=p(a i n i ) + (l-p)(-a i n i ), (4.8) 

where p and 1 — p are the probabilities of being in the two pure states at the two 
ends ±Wj of the diameter whose direction is given by the unit vector n^. 



12 



We can associate an observable a with the real function a(y) = 2aiHi on the 
(n — 1) sphere of points {yi\ykUk = j}- Similarly the mixed state p gives rise to the 
non-negative function P(y) = 1 + 4xiyi with mean value 1. One may then interpret 
P(y) as a sort of classical probability distribution and we have 

(«) = ?/ MvHv)/ [ (4.9.) 

where p is the volume element on the (n — l)-sphere of radius \. If n = 2 the 
probability distribution interpretation is consistent. It is a special case of the general 
quantum mechanical discussion on CP n [14]. 

The difficulty with spinorial Jordan mechanics emphassised by Townsend is 
the introduction of the analogue of a Shrodinger equation. In standard quantum 
mechanics, just as in classical mechanics, every observable generates a flow on the 
space V of pure states which extends to a flow on the space of mixed states M. 
One may alternatively postulate, particularly in the case of open systems, that the 
fundamental equation of motion is a generalized master equation or flow on the 
space of mixed states M. which need not necessarily be induced from a flow on the 
pure states. This is essentally Hawking's viewpoint: he takes M = the space 
of density matrices in standard quantum mechanics. The finite (or infinite ) time 
transformations correspond to his linear $ map called the super-scattering matrix: 

p -> $p. (4.10) 

One says that the $ matrix factorizes if this map may be written as 

p -> SpS^ (4.11) 

for a unitary S'-matrix S C U(n; C). Note that that if one wishes to retain the usual 
interpretation of probabilities relating them to the proportions of the outcomes of 
independent experiments one must have a linear master equation. 

In the case of the standard time evolution for finite dimensional quantum me- 
chanical systems with n states one uses the fact that V = CP n_1 is a symplectic 
manifold and the flow is Hamiltonian. The set of observables generate the group 
U(n; C) which acts transitively on CP n_1 . In the case of spinorial Jordan mechanics 
an observable a certainly gives rise to real valued function (i.e. a^) but V = S n is 
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only symplectic for for n = 2. One way to get a flow would be to take the gradient 
flow with respect to the round metric on S 71-1 . Each such function is a conformal 
Killing potential and so would generate a one parameter subgroup of the conformal 
isometries of S 71-1 . However these do not act linearly on the space of mixed states 
M. = B n . It seems therefore that in spinorial Jordan mechanics one must give up 
the idea that a single observables can generate the dynamics. 

An alternative approach to finding an equation of motion is to use two oberv- 
ables a and b [12]. The analogue of the Heisenberg equations of motion for a density 
matrix p is then taken to be 

p+{a 7/ 9,6} = 0, (4.12) 
where the associator of the three elements {a, p, b} is defined by 

{a,p, b} = (a» p) • b - a« {pm b). (4.13) 

If a, p, b are matrices one has 

{a,p,b}=^[[b,a],p\, (4.14) 

where [6, a] is the ordinary matrix commutator. In ordinary quantum mechanics 
we can always find a and b such that the Hamiltonian H satisfies 

H= l -[a,b] (4.15) 

and so (4.16) reduces to Heisenberg's equation: 

p + i[H,p]=0. (.10) 

A simple calculation shows that (4.12) becomes in components 

Xi = (bidj — ciibj)xj (4-17) 

This is just a rotation in the 2-plane spanned by the n-vectors whose components 
are a\ and bi. It certainly does not correspond to a gradient flow. However if (4.12) 
is generalized to include a sum of associators of pairs of elements one would obtain 
a general element of so(n;R). Thus in spinorial Jordan mechanics the natural time 
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evolution law seems corresponds to rigid rotations of the space V = S n ~ x of pure 
states. 

Following Hawking one might introduce a linear map $ : B n — > B n which 
fixes the totally ignorant state p = -|, i.e which fixes the origin. Such maps would 
be the product of a purity preserving rotation combined with a purity reducing 
contraction. The invertible $ maps then would be the rotations generated as above. 

The conclusion seems to be that in general Jordan's algebraic formulation of 
Quantum Mechanics, in particular the spinorial version, fails the test of providing 
a satisfactory time evolution law unless it coincides with the conventional case. 

Conclusion 

What I have attempted to do in this article is to show that spinorial techniques 
using or Majorana spinors provide a convenient language for describing discrete 
classical systems such as disrete or cellular automata. The observation that spinor 
techniques may be useful is itself is not new. However what is often not greatly 
stressed is that one only needs real numbers. Complex numbers only come in when 
one considers quantum systems and their introduction is not as completely trivial as 
is sometimes supposed. In conventional quantum mechanics the complex numbers 
are needed when considering time evolution. This problem of introducung a time 
function and a complex structure into quantum mechanics becomes especially acute 
in quantum cosmology [15,16] where the consensus appears to be moving towards 
the idea that both are approximate concepts. 

It is therefore of interest to see what form this problem takes in generalizations 
of conventional quantum mechanics. An illuminating example is provided by Jor- 
dan's spinorial algebraic version of quantum mechanics. We found that although 
there is a striking connection between the light cone of n + 1-dimensional Minkowski 
spacetime and spinors. the enterprize founders on the problem of time evolution 
unless it coincides with the orthodox spin | quantum case. However, as I have 
argued elsewhere [15,16], this does not necessarily mean that the complex numbers 
of the spin group Spin(3, 1) = SX(2, C) are to be identified with the complex num- 
bers of quantum mechanics. Indeed one could formulate all of classical physics in 
terms of Majorana spinors and stll have to introduce complex numbers to pass to 
the quantum theory [15,17]. Whilst thay are undoubtedly useful, Weyl spinors are 
in no way obligatory. 
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